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ABSTRACT 



Context. There is observational evidence of propagating kink waves driven by photospheric motions. These disturbances, interpreted 
as kink magnetohydrodynamic (MHD) waves are attenuated as they propagate upwards in the solar corona. 
Aims. To show that resonant absorption provides a simple explanation to the spatial damping of these waves. 

Methods. Kink MHD waves are studied using a cylindrical model of solar magnetic flux tubes which includes a non-uniform layer at 
the tube boundary. Assuming that the frequency is real and the longitudinal wavenumber complex, the damping length and damping 
per wavelength produced by resonant absorption are analytically calculated in the thin tube (TT) approximation, valid for coronal 
waves. This assumption is relaxed in the case of chromospheric tube waves and filament thread waves. 

Results. The damping length of propagating kink waves due resonant absorption is a monotonically decreasing function of frequency. 
For kink waves with low frequencies the damping length is exactly inversely proportional to frequency and we denote this as the 
TGV relation. When moving to high frequencies the TGV relation continues to be an exceptionally good approximation of the actual 
dependency of the damping length on frequency. This dependency means that resonant absorption is selective as it favours low 
frequency waves and can efficiently remove high frequency waves from a broad band spectrum of kink waves. The efficiency of the 
damping due to resonant absorption depends on the properties of the equilibrium model, in particular on the width of the non-uniform 
layer and the steepness of the variation of the local Alfven speed. 

Conclusions. Resonant absorption is an effective mechanism for the spatial damping of propagating kink waves. It is selective as 
the damping length is inversely proportional to frequency so that the damping becomes more severe with increasing frequency. This 
means that radial inhomogeneity can cause solar waveguides to be a natural low-pass filter for broadband disturbances. Hence kink 
wave trains travelling along, e.g., coronal loops, will have a greater proportion of the high frequency components dissipated lower 
down in the atmosphere. This could have important consequences with respect to the spatial distribution of wave heating in the solar 
atmosphere. 

Key words. Magnetohydrodynamics (MHD) — Waves — Magnetic fields — Sun: atmosphere — Sun: oscillations 



1. Introduction 

The first observations of post-flare transversal coronal loop 
oscillations by t he Transition Regi o n and Coronal Explorer 



ing, e.g., footpo int leakage ("B erghmans & de Bruvne', 'l995t 
De Pontieu et al.L [200 1 ). phase mixing (Hevvaer ts & Priest , 



(TRACE) (e.g.. lAschwanden et all 119991; iNakariakov eTal 



1983L iRobertsL l2()02t lO fman & Aschwanden. ' l2002h and res; 
onant absorption ('Ruderman & Roberts, 2002; Goossens et al 



1 1 9991; lAschwanden et al. , 2002), inspired much development in 



I magnetohydrodynamic (MHD) wave theory. This observational 
breakthrough was important since estimated wave parameters, 
such as frequency and amplitude allowed us to implement mag- 
, netoseismological techniques to probe the plasma fine struc- 
I ture of the Sun's atmosphere, an idea initially proposed by e.g., 
iUchida ( 1970) and Roberts et al. ( 1984). It is now commonly ac- 
cepted that these transversal waves are t he kink mode from MHD 
wave theory (see e.g. . .Edwin & Robertslll983|). ahighlv magnet- 
ically dominated, i.e., Alfvenic wave (see lGoossens et aUl2009l 
for a discussion on the nature of kink waves). The observed 
post-flare kink waves in coronal loops have two main defining 
characteristics, firstly they are standing modes and secondly, 
they are strongly damped oscilla tions (in about 1-4 periods, 
see e.g., lAschwanden et all l2003l) . Initially there were several 
physical mechanisms proposed to explain the observed damp- 



Send offprint requests 
j aume . terradas@uib . es 



Terradas, 



e-mail: 



2002) and more recently loop cooling (Morton & Erdelvi, 200' 
Thus far, resonant absorption, caused by plasma inhomogeneity 
in the direction transverse to the magnetic field, has proved the 
most likely candidate in explain ing the observed short damp- 
ing times in coronal loops (see iGoosseni. l2008l for review). 
Consistent seismological studies based on resonant absorption 
using the observed values of periods an d damping times of 
sta nding kink waves were carried out by lArregui et alJ (l2007l) 
and lGoossens et all (l2008l) . These two studies show that, at least 
for a collection of 1 1 loops, resonant absorption provides an ex- 
planation of the observed periods and damping times. The ob- 
servational signatures of the alternative cooling loop damping 
mechanism proposed by iMorton & Erdelvil (l2009i) . differ from 
resonant absorption by the fact that frequency changes as a 
function of time (the resonant damping theory developed so far 
is restricted to static equilibria and therefore frequency is ex- 
pected to be constant in time). Extensive MHD modelling has 
also shown that resonant absorption is the most likely explana- 
tion for the damping of transverse oscillations in prominence 
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fine structures (see 



lArregui & Balleste^ 



Arregui et all I20Q8I: ISoler et all l2009allbl: 
2010h . 

coronal loop kink oscillations are being at- 
process of resonant absorption, one can 



If indeed, 
tenuated by the 
exploit this to estimate the transverse plasma inhomogeneity 
length scales us i ng ob served frequencies and damping times 
dGoo ssens et al.. 2006). The original equilibrium models by 
e.g., Rud erma n & Roberts ( 2002), to study resonant absorption 
consisted of monolithic loop structures. However, there has 
been some observational evidence by e.g.. lAschwandenl (l2005h 
that coronal loops are composed of many different strands, 
possibly at different temperatures and densities. To model 
this loop multi-thread str ucture, there was an initial numerical 
study by iTerradas et al.l (jfeOOS) and it was found that the 
process of resonant absorption was still an efficient damping 
mechanism in more complex and realistically structured loop 
models. Further study into the properties of standing kink 
waves was also undertaken relating plasma inhomogeneity 
in th e direction of the magnetic fiel d caused by e.g., den- 



sity dPfa z et al 



2005t 



'2004' 



Andries et al.. 



2005b; Arregui et al.' 
Pvmova & Ruderman, 2006; Erdelvi & Verth, 2007; 
Verthetal., 2007) or magnetic (IVerth & Erdelvil [20081: 
Ruderman et al., ,2008 ) stratification. It was shown that 
eigenfrequencies and eigenfunctions of coronal loops with 
longitudinal stratification, were altered in such a way that 
one could determine, e.g., the coronal density scale height 
by estimating the ratio of the fundamental mod e to that of 
higher overtones, a technique first proposed by lAndries et alj 
( l2005ah and later developed bv IVerth et aTl ( l2008h to correct for 
magnetic stratification. 

All this theoretical development to describe kink waves in 
coronal loops with realistic plasma stratification in the trans- 
verse and longitudinal direction was restricted to the study of 
standing waves, since this was what was detected in TRACE 
data. However, recently it has come to light that there are 
also ubiquitous small amplitude propagating transversal MHD 
waves in the solar atmosphere. These were first observed by 
the novel Coronal Multi-Channel Polarimeter (CoMP) instru- 
ment (Tomczvk et al., 2007). Moreover, Tomczvk & Mclnto^ 
(|2009|) have been able to separate outward and inward propa- 
gating wave power. It was found that the outward power was 
greater than the inward power by about a factor of two and this 
can only be explained if the waves are damped in situ (see also 
^ascoe et al., 2010). The reason we could not detect these prop- 
agating waves previously with TRACE is that the amplitudes are 
of the o rder 50 km, wh i le the TRACE resolution is only about 
800 km. Tomczvk et all (l2007h originally interpreted thes e wave 
as Alfven waves, however. Van Doorsselaere et a n (l2008h subse- 
quently argued that the observed waves were actually more con- 
sistent with the propagating kink mode. Although both modes 
are dominated by the restoring force of magnetic tension, in the 
geometry of a solar flux tube, e.g., a magnetic cylinder, a pure 
Alfven wave is strictly torsional with no transverse component 
and therefore completely incompressible. On the other hand, a 
kink wave propagating in a flux tube has a transverse perturba- 
tion component and is weakly compressible, at least in the linear 
regime (see e.g., Goossens et al., 2009). 

In the present paper we shall restrict the study to inves- 
tigating the effect of transverse plasma inhomogeneity on the 
propagating kink mode. Although the two problems of standing 
and propagating kink waves are closely related, since a stand- 
ing wave is a superposition of two propagating waves, there are 
some differences that need to be considered. The standing trans- 
verse oscillation is the result of an initial value problem, i.e., an 



initial disturbance in the solar corona, e.g., a CME or a flare, in- 
duces the oscillations of the loops at their natural frequencies or 
eigenmodes. Alternately, the transverse travelling waves have a 
forced nature since the photosphere acts as a driver. The frequen- 
cies of the kink waves observed by CoMP show a peak around 
5 min, indicating a p-mode driven photospheric origin. The spa- 
tial scale of the driver at the base of the tube is also important 
in order to excite kink oscillations. From the properties of MHD 
waves in a flux tube we know that a purely incompressible exci- 
tation excites purely incompressible Alfven waves if the driver 
is strictly localised inside or outside the tube (assuming a homo- 
geneous loop model). On the other hand, an excitation located 
both inside and outside the tube invariably excites transverse os- 
cillations since an almost incompressible surface wave between 
the two media will be established, i.e., the kink mode. 



2. Waveguide model 

To understand the effect of radial inhomogeneity on propagating 
kink waves we consider the relatively simple equilibrium model 
of a cylindrical axi-symmetric flux tube of radius R with a con- 
stant axial magnetic field B, and with a density contrast of pjpe- 
The subindex "i" and "e" refer to the internal and external part 
of the tube, respectively. It is also assumed that the tube has a 
smooth variation of density across the waveguide boundary (lo- 
cated atr - R) with a characteristic spatial scale /. For simplicity 
a sinusoidal density profile connecting the internal and external 
part of the tube is imp lement ed (see e.g., iR uderman & Roberts!, 
I2OO2I; iGoossens et aU. 120021; ITerradas et al.. ,2006) . The role of 
the inhomogeneity at the loop boundary is crucial since this is 
where the process of resonant absorption invariably takes place. 
This basically means that the transverse displacement of the 
whole tube is converted into azimuthal motions localised at the 
tube boundary. Due to this energy conversion, the transverse mo- 
tion is attenuated, while at the same time small scales are cre- 
ated in the nonuniform layer due to phase mixing. Phase mix- 
ing causes a cascade of energy to smaller length scales, where 
the dissipation b e comes more efficient. The re a der is referred to 
iGoossensI (lloosh ; iRuderman & Erdelvil (l2009t) ; iTerradasI (l2009h 
and references therein for further det ails about th i s robus t damp- 
ing mechanism. Recent studies by Soler et aTl (12009 a) in the 
context of modelling kink oscillations observed in solar promi- 
nences (complex coronal magnetic structures with relatively 
cool and dense plasma), show that the process of resonant ab- 
sorption still survives even when the plasma is partially ionised. 



3. Spatial damping in the TT approximation 

Before we embark on the analysis of MHD waves in the thin 
tube (TT) approximation, let us explain what this approxima- 
tion actually is. In what follows we shall consider both standing 
and propagating waves. A standing wave is the superposition of 
two propagating waves travelling with the same frequency and 
wavenumber but in opposite directions. For the standing wave 
problem the axial wavelength (or the axial wavenumber k,) is 
specified and the corresponding frequency is determined. A TT 
in this case means that the axial wavelength is much longer than 
the radius of the tube so that kM 1. Hence the TT approx- 
imation for standing waves is the long wavelength approxima- 
tion. For propagating waves the frequency co is specified and the 
corresponding wavelength is determined. In this situation, a TT 
means that during one period, defined by the frequency oj, a sig- 
nal travelling at the Alfven speed can cross the waveguide in 
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the radial direction many times or loKva/R) ^ 1- Hence for 
propagating waves, the TT approximation is the low frequency 
approximation. The benefit of the TT approximation is that it 
enables us to obtain simple mathematical expressions which are 
very accurate, allowing us to gain a physical insight into the 
problem. 

3. 1 . Homogeneous magnetic cylinder 

For a homogeneous magnetic cylinder (i.e., no inhomogeneous 
layer, / = 0) the dispersion relation for the kink mode (m - I) 
and the fluting modes (m > 1) in the TT approximation is (see 
iGoossens et"Sll2009l) 



Pi(w - 
where 



to' 



2,, 2 



0, 



(1) 



A ^^^BiK^p). (2) 

The dispersion relation Eq. ([T]i specifies a functional dependence 
on frequency u and wavenumber k~. This can be studied for the 
cases of either standing or propagating waves. For the propagat- 
ing wave study, the wavenumber k- is specified and the disper- 
sion relation is solved for frequency o), leading to the well known 
result 

2 Pit^ii+Pet^ie _ 2 
CD = = W,, (3) 

Pi +Pe 

where w» is real. For equal magnetic field strength inside and 
outside the cylinder this expression can be further simplified to 



2Bt 



-K = 



2pi 



2 ; 2 _ 2 



(4) 



P(Pi +Pe) ' Pi +Pe 

For propagating waves we consider waves generated at a given 
location with a real frequency w* and solve the dispersion rela- 
tion for k,, resulting in 



,2 _ /^(Pi +Pe) 2 _ Pi +Pe _ ,2 



V 



Ai 



where K is also real. The solution given by Eq. Q corresponds 
to the well known undamped kink wave (and the fluting modes). 
To derive Eq. ([T]) it is implicitly assumed that ojR/va <k 1 but 
since not all frequencies will satisfy these conditions later we 
present results for any frequency of driver without this restric- 
tion. 

3.2. Inhomogeneous magnetic cylinder 

The thin boundary (TB) approximation means that the non- 
uniformity is confined to [R - 1/2, R + 1/2], with l/R <s: 1, so 
that the non-uniform layer coincides with the dissipative layer. 
This approximation results in the mathematical simplification 
that MHD waves can be described by solutions for uniform 
plasmas that are co nnected over the dissipative layer by jump 
conditions (see e.g.,'Sakurai et al.','l99l'; Goossens et al.,'1995'; 
[rirry & Goossens, 1996). The use of the TB approximation has 
bee n applied in e.g . . iGoossens et al.l ( l2006l) . iGoossensI (l2068h 
and lGoossens et all (l2009h . The inclusion of the effect of an in- 
homogeneous layer is reasonably simple in the case of the thin 
boundary approximation and results in the following dispersion 
relation, 

Pi{a/ - W^i) + Pe(w^ ■ 



in 



Pirh) |Aa 



■pi [u? - w^i) pe (w^ -(A^)- (6) 



Here ta denotes the position of the Alfven resonance. In the 
TB approximation it is natural to adopt - R since l/R «; 1 
and ta &]R - l/2,R + l/2[. Use of the jump condition is not re- 
stricted to the thin non-u n iform layer a s can be see n from e.g., 
iTirrv & Goosseiisl ([19961) : iTirrvetalJ ([19971 [19981) . However, 
this condition requires numerical integration of the ideal MHD 
equations in a non-uniform plasma up to the dissipative layer 
In the present paper we do not intend to use numerical integra- 
tion of the ideal MHD equations relating to thick boundaries. 
However, we shall use the results obtained with the TB approxi- 
mation for thick boundaries for a comparison with the results of 
a full numerical calculation in Section|5] 

Note that the effect of resonance is contained in the right 
hand side of Eq. ((6]). Again we can view the dispersion relation 
Eq. ((6]l as a relation for either standing or propagating waves. In 
the case of standing waves the wavenumber is real {k- - K) and 
the frequenc y is complex. This case has been considered previ- 
ously by e.g. . IGoossens et all d 19921 12009) . Let us now focus on 
propagating waves with a given real frequency (o) - cj,) and 
complex wavenumber The imaginary part of the wavenumber 
indicates that the wave, as it propagates, is damped due to reso- 
nant absorption. Now we assume that 



k^ — kji + iki. 



(7) 



The purely imaginary term in Eq. ((Tj) reflects the damping im- 
posed on the wave and since the damping is in the spatial do- 
main the wavenumber is now complex. If we approximate k^ by 
k^ + 2ik^ki (we assume weak damping, i.e., ki «i ^r) we have 
that A;r X! A;, (given by Eq. ||5l) and after some algebra we find 



b. - ^ ^ ^' ~ P""^ 2 

k, ~ 8/?p(rA)|AA| (Pi+Pe) 

By the fact that, 

,2 



(5) p('"a) |AaI = cj 



we finally obtain the following expression 

ki n m(pi — pe)^ 1 

k^ S R pi+pe \dp/dr\ ' 



(8) 



(9) 



(10) 



Because k- is complex we define the damping length as Ld = 
l/ki while the wavelength is simply A - 2n/K. A useful quantity 
is the the damping per wavelength which is 



A R Pi + Pe 



For a sinusoidal density profile it can be shown that 



-PL 

2 



■Pe 



and Eq. ( ITTT i reduces to 

Lp _ 2 I R pi+pe 
A n m I Pi — Pe 



(11) 



(12) 



(13) 



Equation ([T3T l clearly shows the dependence of the damping per 
wavelength on the thickness of the layer. The wider the layer, 
the stronger the spatial attenuation of the wave. This is not sur- 
prising, since we can relate this result to the expression for the 
temporal damping, i.e., wavenumber assumed real (k, = k,) and 
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frequency complex. The well known formula in the case of the 
temporal damping for a standing wave is 



P n m I Pi - pe 



(14) 



If we compare this expression with Eq. ( fT3l l we note that the 
damping per wavelength for propagating waves and the damping 
per period for standing waves are exactly the same. Hence, in the 
TT approximation spatial and temporal damping are completely 
equivalent. From Eqs. ( flj] ) and (fT4l i we obtain the simple result 



Lb 



P ■ 



(15) 



Propagating kink wav es we re recently studied by 
IVasheghani Farahani et al.l (l2009h in the long wavelength 
limit in X-ray jets in the solar atmosphere. These authors com- 
puted the ratio of the damping time to the period for standing 
waves and used this quantity to discuss the spatial damping of 
propagating waves. The TT approximation result presented in 
Eq. (fTSl l. therefore validates the use of the damping expression 
derived for kink standing waves by Vasheghani Farahani et al.l 
(l2009i) to interpret the attenuation of propagating kink waves. 
The TT damping relations given by Eqs. ( fT3] ) and ( fT4l i have 
further significant consequences. For standing waves we rewrite 
Eq. (O as 



P 
where 



m ' 



ft 



2Rpi+Pe 



TT I pi-pe 



(16) 



(17) 



which only depends on the parameters of the equilibrium model, 
not on the particular type of MHD wave mode defined by the 
value of m (note that for the kink mode m - 1). The period is 
defined by 



27T 
P = —, 



(18) 



and for standing waves oj is related to the wavenumber k, by the 
dispersion relation given by Eq. (|4|, i.e.. 



TT 2pi 

w = n VAi — 



'I A/a 



(19) 



where vai is the internal Alfven speed, n = 1, 2, 3, ... is the lon- 
gitudinal mode number and L is the total length of the waveg- 
uide, e.g., a coronal loop (we have used that k, - mtt/L). 
Equation ( fTSI l can then be written as 



Td 

TAi 



2fe j pi + Pe 1 

m y 2pi n 



(20) 



where TAi = ^/i^'Ai is the Alfven transit time in the longitudinal 
direction. Equation ( l20b has interesting consequences. Firstly, 
it indicates that fluting modes (m > 1) have shorter damp- 
ing times than the kink mode (m - 1). Secondly, it shows 
that the damping time for a standing wave is inversely pro- 
portional to the longitudinal mode number n, i.e., the damping 
time is inversely proportional to the wavelength of the stand- 
ing wave, so that higher overtones (with shorter periods) are 
damped faster than low order overtones, e.g, the fundamental 



mode. Fortunately, there have been some signatures of over- 
tones in coronal loop standing kink waves detected in TRACE 
data (see e.g.lVe rwichte et al., 2004; De Moortel & Bradv, 20071 
IVerth et al.L 120081: IVan Doorsselaere et all 2009). For a particu- 
larly clear example of the first overtone damping before the fun- 
damental mode, that could be explained by resonant absorption 
attenuating the higher harmonic faster, see the Morlet wavelet 
transform in Figure 5 of Verwichte e t al., (.2004.) . 

Now considering the spatial damping of propagating waves 
we can also write Eq. ( fTlT l as 



A 



m 



(21) 



The wavelength is defined as A - 2TT/k, and is related to the 
frequency by the dispersion relation in Eq. (|5]l. Equation (|2TI) 
can then be rewritten in terms of the flux tube radius as 



R 



2tt 



VAi ^£ 



2pi 



(jjR m y Pi -H Pe 



(22) 



where ojRIvpji is a dimensionless frequency (note that the TT 
approximation means a>R/vAi 1). Denoting this frequency as 



ioR 
J ~ ' 

VAi 

then 



2pi 1 



Ld ^ 
R m y^pi+pf, f 



(23) 



(24) 



We shall refer to the relation between the damping length and 
frequency defined in Eq. (l24l i as the TGV relation for propa- 
gating kink waves. In what follows, it will become clear that 
the TGV relation is actually a very good approximation of the 
dependency of damping length on frequency for all relevant fre- 
quencies as illustrated in Fig. 2. The expression given in Eq. ( l24l l 
has important consequences, as it shows that L^jR is inversely 
proportional to / for propagating waves, i.e., the damping length 
is inversely proportional to frequency, so that high frequency 
waves are damped in shorter spatial scales than their lower fre- 
quency counterparts. This means that for driven waves propagat- 
ing upwards from the photosphere each frequency has a different 
penetration height into the solar corona. Therefore, resonant ab- 
sorption provides a natural filtering mechanism for a broadband 
disturbances, e.g., like those observed by T omczyk & Mclntoshl 
(l2009h . with lower frequency waves being least affected by the 
damping process, propagating to higher heights in the solar 
corona. 



4. Spatial damping beyond the TT approximation 

The TT approximation is useful because it makes the dispersion 
relation analytically solvable. In Table [T] we list observed esti- 
mates of kji for various MHD wave modes and it can be seen 
that for many solar atmospheric waves, e.g., standing kink waves 
in coronal loops observed by Aschwanden et al. ( 2002), the TT 
approximation is reasonably valid . However, for kink waves in 
filament threads dLin et al 1 120091) or torsional Alfven waves in 
chromospheric magnetic bright points (iJess et al 1 I2OO9I) may 
have values of kji ^ 1. To address the kink waves observed 
in this regime, it is relatively straight forward to relax the TT ap- 
proximation and to take into account the effect of the finite tube 
radius on damping. 
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Table 1. Estimated range of A: in observed propagating (P) and standing (S) MHD wave modes. 



Reference 


MHD wave mode 

itifpmrptatioTi nf aiitfior^ 


Width of 

wjivpcTiii Hp 

W £1 V eg Ul \J\^ 

(Mm) 


Wavelength 
(Mm) 


LR 


1 


Kink (S) 


5.9- 11.5 


328 - 552 


0.03-0.11 


2 


Alfven (P) 


0.2 


> 4.0 


< 0.16 


3 


Kink or sausage waves (S) 


0.4 


2.5-3.1 


0.41-0.5 


4 


Torsional Alfven (P) 


2.0 


4.5 - 6.0 


1.04-1.39 


5 


Kink(P) 


0.2- 1.0 


3.4-12.9 


0.03 - 0.9 


6 


Kink(P) 


0.43 - 0.66 


> 250 


< 0.12 


7 


Alfven (P) 


9.0 


> 393 


< 0.07 


8 


Kink or Alfven waves (P) 


9.0 


> 180 


< 0.16 



Re ferences. (l)lAschwanden et al. (2002); (2) De Pontieu et al. (2007); (3) Fuiimura & Tsunetj ( 120091) ; (4) lJess et alj ( 120091) ; (5) lLin et al.l ( 120091) ; 
(6) IOkamoto et alJ (l2007h ; (7) rTomczvk et al. (.2007) ; (8) .Tomczvk & Mcintosh 1.2009lv 



4. 1 . Homogeneous magnetic cylinder 

For a homogeneous magnetic cylinder the dispersion relation, 
taking into account a finite radius, is 

Piico^ - col^) + p,((J - cjI,)F(co, k,) = 0, (25) 



where 



Fico, L) = - 



k J'JkR)K„{k,R) 
k, J„(kiR)K;„ik,Ry 



- 4i ,2 

2 ' e 



(26) 



(27) 



Note that the function F depends both on frequency oj and wave 
number k^, taking into account that the tube has a finite width. In 
the limit of a TT, F = I (the Bessel functions are approximated 
by their small arguments expressions) and we recover the disper- 
sion relation given by Eq. ([T]i. For a standing wave fc; is fixed and 
Eq. ( l25b is solved for co. For a propagating wave cj is fixed and 
Eq. ( |25] ) is solved for k^. The solution to Eq. ( |25] ) corresponds 
again to the undamped kink wave. To solve this equation ana- 
lytically for a finite width tube one would have to use additional 
terms in the asymptotic expansions of the Bessel functions, com- 
plicating matters. To avoid these cumbersome calculations we 
simply solve Eq. ( 1251 ) numerically, with the solution differing 
from the solution given by Eq. (|5]l outwith the TT regime. 

4.2. Inhomogeneous magnetic cylinder 

When we add a thin non-uniform layer we obtain the complex 
dispersion relation 

Pi(a/ - aj\-) + pja? - a)\j F(a), k.) = 



in- 



m/R 



p(r a)\Aa 



■Pi (w^ _ ^2 ^^^^ ^^2 _ ^2^^ ^(^^ 1^^-^^ (28) 



where 



G(w, k^) = - 



m 1 K„,{k,R) 
Rk,K'Jk,Ry 



For the TT we can again use the asymptotic expansions of the 
Bessel function so that G s 1, hence in the TT limit, Eq. (|28]) 
is simplified to Eq. (|6]l. Equation ( 1281) can be solved for a com- 
plex frequency u) with a specified real wavenumber K or for a 
complex wavenumber k, for a specified real frequency w*. By 
solving Eq. ( l28b in the case of a propagating wave with real fre- 
quency we find that the real part of complex h is K (we get 
again Eq. ||25]| . note that now K is different from Eq. ||6l), while 
the imaginary part is given by 



T 



(30) 



where 



n m 

r = - 



(Pi - Pe) 



^R\pi+ p,Ficj„ k,)] [1 + F{aj„ k,)f 
\dp/dr\r^ 



and 



= 1 + 



Pi -Pe dF_ 

2 [pi + Pe/^(w„ k,)] [1 -I- F(w., k.)] dk 



(31) 



(32) 



The expressions are slightly more complicated than in the TT 
approximation (see Eq. IfTOll ). but once k, is known the different 
terms can be easily calculated. The corresponding damping per 
wavelength is 



A 



(33) 



In the TT limit, i.e., when F{(jl>^,K) and G(a)t,kt) 1 and 
dF/dk —> 0, we recover Eq. ( fTTT ). For a sinusoidal density profile 
we simply have to make use of Eq. (fT2l) in Eq. (|3TI) . Now let us 
define 



(29) = 1 



CO, 



Pi 



dF 



2 \pi+p^F{a)„K)\[l+F{a).„K)] dco 



(34) 
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As in the previous Section we compare the damping per wave- 
length for propagating waves with the damping per period, 
which in the non-TT approximation is given by 



12. 
P 



]_N 
2^Y' 



(35) 



By Eqs. ( l33T l and ( l35T l. there is a clear correspondence between 
the temporal and the spatial damping beyond the TT approxima- 
tion. Furthermore, a general expression that relates the temporal 
damping of standing waves and the spatial damping of propa- 
gati ng waves can be deri ved following the analysis of Appendix 
A in lTagger et alJ (1 19951) . We can write the complex dispersion 
relation given by Eq. ( l28l l as 



Dr(w, k) + iDiioj, k) = 0. 



(36) 



To make it clear, in Eq. (1361 ), Dr is equal to the LHS of Eq. ( l28T l 
and Di equal to minus the RHS of the equation. In the case of 
spatial damping, if <y» and k^ are the solutions of Z)r(w, k) - Q 
then it is easy to see by making a Taylor expansion around the 
solution that 



D, 



dD^/dk 
while for temporal damping 
Di 



Lull 



Combining these two expressions we find that 



A 



h dD^/dk td 



w. dD^ldco P ' 
which is simply (see also lPascoe et al 



A 



Vph P ' 



(37) 



(38) 



(39) 



(40) 



This result is valid when the damping is not too strong (see e.g., 
iTag^er et al.L Il995h . an assumption that we have akeady made 
in the derivation of the damping per wavelength {ki ^r). In 
the TT approximation kink waves are weakly dispersive so that 
Vgi a: Vph hence, in agreement with the results found in Section|3] 
the damping per wavelength is exactly the same as the damping 
per period but in the regime when the TT approximation is not 
applicable the group speed can differ from the phase speed. 

From the previous results we can also derive alternative for- 
mulae for the damping per wavelength and damping per period 
in terms of the imaginary part of the dispersion relation. 



A 



1 



dF 



dk 



(41) 



12. 
P 



L0» 



1 



In Di(cj^,k^) 



dF_ 

d(j) 



(42) 



Using Eq. ( l40b we can identify in these expressions the terms 
related to the phase and group speed. 
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Fig. 1. Damping per wavelength as a function of the dimensionless fre- 
quency (/ = cjR/vAi) for three different widths of the inhomogeneous 
layers. The solid line corresponds to the analytical results, the dashed 
line represents the full numerical solution of the resistive eigenvalue 
problem. The dotted line corresponds to the TT approximation, valid 
when / ^ 0. In this plot pi/pc = 3. 



Next we study, using the analytical results given in this 
Section, how the damping per wavelength depends on the dimen- 
sionless frequency of the driver/, defined by Eq. ( l23l l. for partic- 
ular cases (see FigurelTJ. We see that the wider the layer the more 
efficient the attenuation (smaller damping per wavelengths), in 
agreement with the analytical results in the TT approximation. 
For small / the damping per wavelength tends to the TT value 
(see dotted line). The value of Lu/A increases monotonically 
with a>. The deviation with respect to the TT results is smaller 
for thick layers. In the TT approximation, which for propagating 
waves, is the low frequency approximation, Ln/A is indepen- 
dent of frequency, since this approximation does not take into 
account the variation of frequency, i.e., the frequency is only 
presumed to be small. The analytical results for Lq/A beyond 
the TT approximation, take into account the dependence on fre- 
quency. The value of Lu/A now undergoes a moderate increase 
when we move from low to high frequencies. However, the re- 
ally interesting quantity to calculate is the damping length itself. 
In order to make the dependence of Lq on frequency more ex- 
plicit we follow the same line of reasoning in Section 13.21 and 
rewrite Eq. (fTsT l as 



(43) 



where the quantity ^ew now depends on the parameters of the 
equilibrium model and the characteristics of the wave itself. 
Again A = 2n/k, and k- is related to the frequency by the dis- 
persion relation Eq. ( |25] |. Since we have abandoned the TT ap- 
proximation, a simple analytical formula that relates k, to oj is 
not readily available but we can always write 



k,R^f<P(f), 



(44) 



where / is the dimensionless frequency defined in Eq. ( |23] ) and 
(/'(/) is a function that we can determine numerically. Hence, 
implementing the function ifr(f), we have that 



(45) 
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which in the low frequency limit is equivalent to the TGV rela- 
tion given by Eq. (l24l ) for m = 1, i.e., 



^EW ^ and ^ ^ y (46) 

as / — > 0, where is defined by Eq. ( fTTl i. In Eq. ( l45T l, govern- 
ing the damping length for propagating waves, since </'(/) and 
^Ew are slowly varying functions of /, the main dependence of 
Ld on / is contained in the factor 1 //, as in the low frequency 
limit shown by Eq. UM . This is illustrated in Figure |2] where 
Lx)jR is plotted as a function of / for several values of IjR. The 
dependence of the curves on \lf \& very clear and the non-TT 
and TT solutions (compare the solid with the dotted lines) tend 
to overlap in the limit of / — > 0, where both cases are accu- 
rately described by the TGV relation. Note that even for / — > 1 
the TGV relation still describes quite well the behaviour of the 
damping length with frequency. Again, the potential of reso- 
nant absorption as a frequency filter is clearly demonstrated in 
Figure |2] with high frequency waves being damped in shorter 
spatial scales than low frequency waves (for fixed ///?). 



1QQ0 r~\ \ • r 




ol . . . \ . . . \ . . . \ . . . \ . . . I 

0.0 0.2 0.4 0.6 0.8 1.0 



/ 

Fig. 2. Damping length normalised to the loop radius as a function of 
the dimensionless frequency (/ = coR/v/^J for three different widths 
of the inhomogeneous layers. The solid line corresponds to the analyti- 
cal results, the dashed line represents the full numerical solution of the 
resistive eigenvalue problem. The dotted line corresponds to the TT ap- 
proximation calculated using Eq. d24t . valid when / — > 0. In this plot 
Pi/Pc = 3. 



5. Resistive calculations 

The results based on TB approximation, described in the previ- 
ous Sections, are comp ared with the full resistive calculations 
fsee lTerradas et al.l 1200 6. for the details about the method). This 
comparison is useful since in the resistive eigenvalue problem 
there are no implicit assumptions about the TT or TB approx- 
imation, i.e., the TT and TB approximation are not used. The 
resistive calculation applies to any frequency whether it is small 
compared to v^i/R or not and also to equilibrium models that 
have a thin non-uniform layer or are fully non-uniform. We have 
solved the eigenvalue problem (for co) and have used Eq. ( l40t to 
translate from temporal damping to spatial damping. Due to the 
inclusion of resistivity, the eigenvalue problem for the wavenum- 
ber is more difficult to solve than the eigenvalue problem for 
the frequency. The results of the resistive calculations are shown 



in Figures [T] and |2] where the damping per wavelength and the 
damping length are plotted (see dashed lines) as a function of 
the frequency of the driver The agreement between the resis- 
tive computations and the analytical or semi-analytical methods 
is very good. The small differences are due to the fact that in 
the analytical approximations we have assumed that the reso- 
nance is always located at r - R and we have calculated the 
derivative of the density at this position. This explains the small 
deviations from the resistive estimations. A more precise deter- 
mination could be done calculating the exact location of the res- 
onance and then using a slightly modified version of Eq. ( l28l l. 
but since the analytical approximations that we have already de- 
rived are quite satisfactory there is no pressing need to explore 
this issue further. 

6. Conclusions and discussion 

The spatial damping due to resonant absorption of driven kink 
waves has been investigated. The main conclusion of the work is 
that the damping length of propagating kink waves due resonant 
absorption is a monotonically decreasing function of frequency. 
The TGV relation for kink waves was derived, demonstrating 
that for low frequencies the damping length is exactly inversely 
proportional to frequency. In the high frequency range the TGV 
relation continues to be an excellent approximation of the actual 
dependency of the damping length on frequency. Certainly, for 
all physically relevant frequencies the dependency of damping 
length on frequency is accurately described by the TGV relation. 
This dependency means that resonant absorption is selective as 
it favours low frequency waves and can efficiently remove high 
frequency waves from a broad band spectrum of kink waves. 
This has great significance for solar atmospheric kink waves, 
since high frequency waves will tend to lose more power than 
their low frequency counterparts before reaching high altitudes 
in the solar corona, with the exact percentage power loss depend- 
ing on the properties of the equilibrium, in particular the width 
of the non-uniform layer and steepness in variation of the local 
Alfven speed. With respect to mode conversion, the process of 
resonant absorption will cause the higher frequency waves to be 
attenuated more because the global kink mode will be converted 
into localised Alfvenic modes at lower heights. If the energy of 
these Alfvenic motions is eventually dissipated, then resonant 
absorption should produce a characteristic distribution of the en- 
ergy as a function of height in the solar atmosphere. This could 
have important consequences with respect to the spatial distribu- 
tion of wave heating in the solar atmosphere. 

It was also shown that spatial and temporal damping are ba- 
sically equivalent. In the TT approximation the damping per pe- 
riod and the damping per wavelength are exactly the same. The 
differences in these two quantities arise in the regime where the 
TT is not valid, but even in this situation it is easy to relate the 
spatial and the temporal damping rates through the group and 
phase speeds of the kink MHD waves. This allows us to trans- 
late the results from the temporally damped waves (w complex, 
k real) to spatially attenuated waves (a> real, k complex) due to 
resonant absorption. This mechanism requires the frequency of 
the driver to be between the internal and the external Alfven fre- 
quency of the tube. This might seem a very restrictive condition, 
but in fact it is just the opposite. In the driven problem the fre- 
quency is fixed but the system chooses the proper wavelength 
(along the waveguide) to accommodate the kink mode in the 
tube. This kink mode generated at the base of the loop propa- 
gates upwards along the tube and at the same time is attenuated 
due to the inhomogeneity at the tube boundary. 
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An interesting result is that both the damping length in the 
spatial problem and the damping time in the temporal problem 
are always smaller than in the TT approximation (when / — > 0), 
meaning that waves with short wavelengths or large frequen- 
cies are always more efficiently damped. Note that the observa- 
tions of standing kink waves observed with TRACE and for the 
propagating kink waves detected with the CoMP instrument are 
precisely in the regime where the TT is applicable, i.e., where 
the waves are less affected by resonant absorption. The fact that 
overtones of standing kink waves and high frequency propagat- 
ing waves have proved difficult to detect may be a direct conse- 
quence of the filtering due to resonant absorption. From a differ- 
ent perspective, we have also shown that the damping per wave- 
length (and the damping per period) has a weak dependence on 
the frequency. 

It is necessary to remark that our results are based on a sim- 
ple magnetic flux tube model, i.e., a straight cylinder, with no 
gravity and pressure. The only inhomogeneity in this model is 
plasma density in the radial direction, however this a character- 
istic property of solar waveguides that is observed at all atmo- 
spheric heights, e.g., chromospheric magnetic bright points or 
coronal loops. Thus the theory of resonant damping of propa- 
gating kink waves due to radial plasma density inhomogeneity 
offers a natural explanation f or the dissipation observed by e.g., 
iTomczvk & Mclntoshl (12009 ). However, a detailed comparison 
between the observations and the expected frequency dependent 
response due to resonant absorption is needed to quantify the 
precise spatial distribution of wave heating due to this mecha- 
nism in the solar plasma, e.g., kink wave dissipation as a function 
of both frequency and height. This problem will be the subject 
of a future work (Verth et al. 2010, in preparation). 
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